Generalized measurement of the non-normal 
two-boson operator Z 7 = a\ + 



Matteo G A Parifi 12 , Giulio LandolfS 3 4 , Giulio SolianiJ 4 



^ipartimento di Fisica dell'Universita di Milano, 1-20133 Milano, Italy 
2 Institute for Scientific Interchange Foundation, 1-10133 Torino, Italy 
3 Dipartimento di Fisica dell'Universita del Salento, 1-73100 Lecce, Italy 
4 Istituto Nazionalc di Fisica Nucleare, Sezione di Lecce, 1-73100 Lecce, Italy 

PACS numbers: 03.65.Ta, 42.50.Xa 

Abstract. We address the generalized measurement of the two-boson operator 
Z 7 = a\ + which, for 7I 2 ^ 1, is not normal and cannot be detected by a 
joint measurement of quadratures on the two bosons. We explicitly construct the 
minimal Naimark extension, which involves a single additional bosonic system, 
and present its decomposition in terms of two-boson linear SU(2) interactions. 
The statistics of the measurement and the added noise are analyzed in details. 
Results are exploited to revisit the Caves-Shapiro concept of generalized phase 
observable based on heterodyne detection. 



The two-boson operator 




is normal [Z 7 , ZU = 1 — \ j\ 2 for I7I = 1. In this case the real Xy = ^{Z 7 + Z^) and 
the imaginary Yy = \{Z 1 — Z^) parts of Z 1 commute [Xy, Y 7 ] — and can be jointly- 
measured. Actually they correspond to the canonical sum- and difference-quadratures 
of the two modes 




where, for k = 1, 2, 




1 



(3) 



On the other hand, for I7I ^ 1, we have 
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where x^^ — -^(a^e^ + a^e is a rotated quadrature of the fc-th boson and 

6 1 = arg7. In this case, the two operators do no commute [Xy, Y y ] = |(1 — I7I 2 ) and 
a generalized measurement should be devised. Indeed, the eigenstates of Z 1 for 7 ^ 1 

|*» 7 =!>(«) ® 1 1 7 » 

where D(z) — exp{za\ — z*ai} is the displacement operator and I7)) = 
\J\ — I7I 2 7™|n) ® \n), do not provide a resolution of the identity, we have 

^|z)> 77 ((z| = (l-|7| 2 )|7| 2ata . 

7T 

We first notice that Z 1 = Rg Z^Rg y where R<p = exp(i <f> a3,a 2 ) & n d therefore, 
without loss of generality, we may restrict attention to the case of real positive 7. In 
this case we have 

1 1 
^7 = ^|(9i+792) Y 1 = — (p! - 7p 2 ) (5) 

In addition, we notice that, up to a permutation of the mode labels, Z^ = jZ _ ± 
and therefore, since the multiplicative constant does not influence the measurement 
scheme, we may further restrict attention to the case < 7 < 1. 

The operator Z 1 is defined on the Hilbert-Fock space of two harmonic 
oscillators. A Naimark extension for the operator Z 1 is a triplet (7i a ,T 7 ,<7), where 
T 7 is an operator defined on an extended Hilbert space H12 <8> 7~ia and a is a state 
(density operator) in Ti a , such that for any state R € 'Yin we have 

Tr i2 [R X 7 ] = Tr 12o [R ® <J Re T 7 ] 

Tr 12 [R Yy] = Tr 12a [R ® a Im T 7 ] . (6) 

Equations (J6]) are usually summarized by saying that the operator T 7 traces the 
operator Z 7 . Of course, Eqs. © do not hold for higher moments: the generalized 
measurement of Z 1 unavoidably introduces some noise of purely quantum origin. In 
general we have 

Tr 12 [R X™] ^ Tr 12a [R®a (Re T 7 ) n ] n > 2 

Tr 12 [i? F 7 "] ^ Tr 12a [i? ® a (Im T 7 )"] n > 2 . (7) 

In this communication we look for a minimal Naimark extension, that is an extension 
involving a single additional bosonic mode 03. In general, for operator of the form 

T 7 = Z y + 7(03,03) the trace condition of Eqs. ^ require Tr a <r 7(03,03) = 0, 

whereas the constraint of normality can be written as 

0= [T 7 ,Tt] = [Z^Z\] + [/(a3,4),/(a 3 ,4) + ] • (8) 

It is straightforwardly seen that / (03,03) = ^03 or 7(03,03) = KOg, where k is a real 
constant, are solutions of Eqs. ([6|) and (|8|). In the following we analyze in details 
whether this kind of extensions can be implemented using only bilinear interactions 
among the three modes followed by measurement of quadratures at the output. 
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The measurement scheme is the following: the modes interact each other via 
the unitary operator £/ 7 , which impose the linear transformation 

Ai \ ( ax \ 

U 1 = M | a 2 | (9) 



and then, at the output, the quadratures 

Qi = ^(A 1 + A{) P 2 = -L(Al-A 2 ) 

are measured with the aim of obtaining, upon the definition T 7 = Qi + iP 2 , 

Tr 12 [RX 7 ] =Tn 2a [R®o-Qi] (10) 
Tr 12 [R Y y ] = Tr 12a [R® a P 2 ] (11) 

for any R, and at least one a such that Tr[ca3] = 0. A suitable evolution operator 
Uy corresponds to the transformation 




M = — 1 -7 -k . (12) 



Upon imposing the constraint of unitarity, i.e [Aj, A\] — 5jk, we have the solution 



k= Vl-7 2 , mi = 0, (13) 

m 2 = -^2(l-7 2 ), m 3 = \/27, (14) 
which makes M a U(3) transformation and leads to 

Qi = A= (?i + 792 + V 1 - 7 2 93) , -Ffe = -ji= (f>i - 7P2 - V 1 - l 2 Pzj , (15) 

and, in turn, to T 7 = ai + 70^ + Kag. Notice that no unitary solution can be found 
(for I7I < 1) for the case f(as, a 3 ) = kcl^, i.e. for linear transformation expressing the 
output modes (A\,A 2 , A 3 ) as a linear combination of (ax, a 2 , a\) 

A question arises on how the unitary C/ 7 can be implemented in practice, as for 
example in a quantum optical setting. As it is well known, any SU(3) transformation 
may be decomposed into a set of SU(2) transformation [5J. In our case the U(3) M- 
transformation may be decomposed using three SU(2) transformations followed by a 
7r-rotation. In Fig. [T]we report the explicit decomposition of M. The circle denotes 
a 7T- rotation on the second mode i.e. a unitary of the form R 2 = exp{iTra\a 2 } ■ 
The boxes correspond to SU(2) rotations i.e. to evolution operators of the form 

Bjk{0jk) — exp |— i8jk (o-ja\. + aka^j |, corresponding to the transformations 
,t m..\( "1 \o..ta..\- I cos9 <J sin8 « 



{ Z - 1 c„ s % ){Z) <16 > 

% Actually, a solution involving a SU(1,1) interaction between and 0,3 followed by a SU(2) 
interaction between a± and 0,2 may be found for I7I > 1 and then extended to the whole range 
of I7I by rescaling. However, this solution unavoidably introduces a larger amount of noise compared 
to that of Eqs. II28I I and II29I I and it will not be considered here. 
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Figure 1. Block diagram of the decomposition of the M transformation of 
Eq. il 1 21) into three SU(2) transformations, each involving two of the modes, 
plus a 7r-rotation. The boxes corresponds to evolution operators of the form 

B jk (0jk) = e- lB ^ a ^ +ak ^) (see text). 



By explicit construction we have 



= [Ii <g> R 2 ® I 3 ] [£23(023) <g I 3 ] [B 13 (9 13 ) (g) I 2 ] [B 12 (6 12 ) ® Ii 



where 



cos #23 = 



I + 7 2 



C0Sfyi3 = 



2 7 2 



1 + 7" 



COS #12 = 



I + 7 2 



(17) 



Other decompositions may be also found, allowing for permutations of modes and 
different rotations. For 7 — > 1 the mode 0,3 decouples from the other two modes and 
the scheme reduces to the joint measurement of quadratures for the normal operator 
Zi®. 

Each outcome from the joint measurement of the quadratures Q\ and P 2 
corresponds to a complex number r = Qi + iP 2 that represents a realization of the 
observable T 7 . The probability density of the outcomes K 1 (t) for a given initial 



preparation R 
function 3(A) 



er is obtained as the Fourier transform of the moment generating 



where 



, d 2 \ 
K i(T) = / -r< 



A*t— At* c 



5(A) 



S(A) = Tr i? <g> ct e 



(18) 



(19) 



Using Eqs. {15]) we have exp{AT,j - A*T 7 } = £>i(A) ® D 2 (-A7) ® D 3 (-\n) where 
Dj(z) is the displacement operator for the mode Oj. Therefore, the moment generating 
function rewrites as 



S 7( A ) = Xi2(A) x 3 (-Ak) , 



(20) 



where Xi2(A) = tr [R £>i(A) ® D 2 {— A 7 )] and X3(- z ) = Tr[er £(3(2:)] is the characteristic 
function of the mode 03. Using (|2"0|) is easy to see that the probability density of the 
outcomes is given by the convolution 



KJr) 



1 



H, ( (t)*W 3 (-t/k), 



(21) 



W 3 {z) being the Wigner function of the mode 03, * the convolution product, and 
Hry[z) the density obtained by the Fourier transform of Xi2(A). In turn, for factorized 
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preparations R — Qi ® g 2 the moment generating function X12W = XiW X2(~^) 
factorizes into the product of the characteristic functions of Qi and g 2 respectively, 
and the density H 1 (t) reduces to the convolution of the Wigner functions of the two 
input signals 



1 

-2 



Wi(t)*W 2 (-t/j) 



(22) 



Using (|15|) it is straightforward to see how the variances of the measured quantities 
Qi and P 2 are related to the variances of the quadratures of interest. We have 



AQ^AX 7 2 + i(l- 7 2 )Aq 3 2 
AP| = Ar 7 2 + i(l- 7 2 )A P 2 , 



(23) 



Tr [a g|] and analogously Ap| = Tr[crp|] (remind that Eq. ^ implies 



where Ag| 

Tr[crg3] = Tr[<rp3] = 0). Notice that the added noise in Eq. (|2"3")1 is the minimum noise 
according to generalized uncertainty relations for joint measurement of non commuting 
observables [TJ [5J [21 El E] . On the other hand, the covariance between the measured 
quadratures i.e. the quantity 



S Ql p 2 = 2 Tr i2a [R ® (QiP 2 + P2Q1)] - Tr 12a [R®a Q x ] Tr 12a [R®a P 2 ] , (24) 



may be written as 



s QiP 2 = ^x,,Y,, ~ -(1 - 7 )Tr Q 



2°" (P393 + 93P3) 



(25) 



where Y" 7 = ^Tr^ [-R (X 7 F 7 + ^ 7 X 7 )] - Tri 2 [R X 7 ] Tr i2 [R Y 7 ] is the covariance 
of the desired quadratures. 

Notice that the added noise to the covariance, Eq. (|25|) . may vanish for some 
preparation of the state a whereas the added noise to the variances, Eq. (|23|) . cannot 
vanish for any physical preparation a. This raises the question of the consequences 
of different field states on the statistics of the measurement and, in turn, of the role 
played by preparations of states in concrete experiments. On the other hand, within 
experimental frameworks, one may take full advantage of possible freedom in preparing 
some of the modes. This is definitively the case of the Naimark mode 03, even though 
its preparation needs to be compatible with the prescription (jGj) for the expectation 
values of position and momentum operators. In particular, a valid Naimark extension 
can be obtained by preparing the mode 03 in the vacuum state a = |0)(0| to let its 
contribution to the noise in formula (|25[) to vanish, since Tr a [cr(<73P3 +^393)] = 0, and 
to minimize Aq| and Ap 2 in (|23|) . since both the terms would be equal to one half. 



Each of the other two fields may be, for instance, in one among the most meaningful 
types of states, such as number states, coherent states, thermal states or phase states 
(i.e. eigenstates of the operator C + iS, where C and S are "cosine" and "sine" 
operators respectively) or prepared in an entangled states. If we consider the fully 
separable state described by the density operator g — R (g) a — gi <g) g 2 ® a , where 
gk, with k = 1,2, denotes the preparation for the k-th bosonic field in the arbitrarily 
mixed state gk = X)m=o^™'l TO )( m l on ^ ne Hilbert space Hk, then the system moment 
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generating function is easily obtained by resorting to 



oo 

I 



Tr k [g k D k (a k )]=e-^- J>«L m (|a fc | 2 ) , (26) 



m=0 



where the £ n 's are Laguerre polynomials. For instance, for coherent and phase states 
Eq. (|26| should be used with 



,2ia |2m 



p « =e -|^^L_ and p « = (l- N 2 )M 2m (27) 

respectively (phase state formulae can be used even when dealing with thermal states 
upon the identification z — exp[— i/3Sw], (3 being the inverse of temperature). Suppose 
no specific conditions do constraint, in principle, the preparation for the mode a^- 
Once again a vacuum choice may be advantageous in some respects. Let us therefore 
focus on the specific case of the measurement of Z 1 on the class of factorized signals 
described by R = gi ® |0)(0| where Q\ is a generic preparation of the mode ai while 
|0) is the ground state of the mode a 2 . In this case g = gi g> |0)(0| ® |0)(0|, Eq. 
(|21jl becomes a Gaussian convolution and the moment generating function becomes 
independent of the parameter 7 

S(A)= X i(A)cxpf-i|A| 2 ') . (28) 



2 



The measured variances are thus given by 



AQ 2 = i(Ag 2 + l) AP 2 = i(A P 2 + l) (29) 

Equations (|28[) and (|29[) contain a remarkable result that may be expressed as follows. 
The measurement of Z 1 on the class of states R—g\® |0)(0| does not lead to added 
noise with respect to the measurement of the normal operator Z\ . 
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Figure 2. The scheme for heterodyne detection. 



This result finds a natural application in the context of heterodyne detection, 
where currents of the form (TTJ) show up. As it is known, in heterodyne detection 
a single-mode signal field E\ of nominal frequency u>\ is mixed through a beam- 
splitter with a local oscillator field El whose frequency u>l is slightly offset by an 
amount u)j <C oj\ from that of the input signal, i.e. u>i = u>l + ^>i- A photodetector 
is placed right after the beam-splitter (see Fig. [2]). The output photocurrent, 
which generally depends on fields parameters and on specific assumptions on the 
apparatus, is filtered at the intermediate frequency u>i. In standard optical heterodyne 
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detection (see e.g. |3|), measuring the filtered photo-current corresponds to realize 
the quantum measurement of the normal operator y = ai + a\ [8J, where a\ (res. 
a\) denotes the photon annihilator (resp. creation) operator for the input (resp. 
image) signal. Measuring the real and imaginary parts of the (actually rescaled) 
output photocurrent thus provides the simultaneous measurement of both input field 
quadratures. Nevertheless, it has been also argued that whenever one is not restricted 
to an input field frequency in the optical regime, but, rather, one is concerned with 
microwave (or radio) heterodyning, then the interaction of the input signal field with 
the apparatus of Fig. 2. (approximatively) results in the measurement operator yc = 

yftl + %)*l + W) a l ( see M and discussion in @). Since [y c ,y c ] = 2^ ^ 0, 

Caves measurement operator yc is not compatible with simultaneous measurements 
of signal quadratures. In other words, standard heterodyne detection cannot achieve 
the measurement of the Caves operator and a question arises on whether simultaneous 
phase and amplitude measurements may be accomplished in this case. The answer 
may be found in the results reported above. In fact, the measurement of the Caves 
operator corresponds to the generalized measurement of the non-normal operator 



Z lc = ai + 7c 4 . 7c = J LUl , UJI < 1 (30) 

In the light of our previous results, we thus learn that the simultaneous measurement 
of the field quadratures for a quasi-monochromatic signal can be realized even in the 
case when the heterodyne apparatus yields a measurement operator of the Caves type, 
Eq. (J3DJ) - T° this aim, it suffices to generalize the heterodyne detection scheme by 
introducing a single boson Naimark mode and letting it interact with the other modes 
through the linear transformation Moreover, a suitable preparation enables one 
to avoid additional noise with respect to that resulting in the measurement of signal 
field quadratures within the framework of the standard optical heterodyne detection. 

It is worth also discussing the matter from the point of view of phase operators 
since our results can be used to proceed in defining a feasible phase within the Caves 
description of heterodyning. Since the operator T is normal, then its associated self- 
adjoint phase operator 

1 T 

6> T = _ In— (31) 

can be defined unambiguously indeed so that cosine and sine quadrature operators 
c = 1 ^e T + e -ie T j 5 s = J_ ^e T _ 

obey the correct relation C 2 + S 2 = 1. It is now in order to recalling that the 
two-modes relative number state representation discussed by Ban (see [TT] and Refs. 
therein) fits fairly with the feasible phase concept of Shapiro and Wagner (namely, 
the shift phase operator associated with the Shapiro- Wagner measurement operator 
y = a± + c4). Upon defining the 3-mode relative number operator N = Ni — N2 — N3 , 
where Nk — o^Ofc (k = 1, 2, 3), one gets 

[e ldT , N] = e if>T , [N,6 T ]=i. (32) 



These relations are what one expects for genuine phase operators. In other words, a 
feasible phase can be naturally defined even in the Caves description of heterodyning at 
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the cost of introducing of a Naimark mode and generalizing the 2-modes relative state 
representation to a 3-modes one. The commutator [N, 9t) can then be interpreted 
as the canonical conjugation of the feasible phase for Caves heterodyne measurement 
operator with respect to the operator mode number difference N. 

As final comments, notice that tracing out the Naimark mode a^, and introducing 
symmetric ordering when needed in Eqs. (|3ip - (|32p . formulae given in |10j are 
recovered. Further, it would be of interest to move towards the direction of generalizing 
the relative number state representation for the description of the phase operator of the 
generalized heterodyne measurement we have introduced in this communication, and 
more generally for operators describing linear amplifiers involving more than three 
modes. This is also concerned with the investigation of the possibility to extract 
basic algebraic structures underlying these systems to generalize algebras given in [10] . 
These issues are currently under investigation and results will be reported elsewhere. 

This work has been supported by MIUR through the projects PRIN-2005024254-002 
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